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Abstract 

Using the NMC and E665 nuclear structure function ratios F^/F® and F^/F^ from deep in- 
elastic lepton-nucleus collisions, and the E772 Drell-Yan dilepton cross sections from proton- 
nucleus collisions, and incorporating baryon number and momentum sum rules, we determine 
nuclear parton distributions at an initial scale Q$. With these distributions, we study QCD 
scale evolution of nuclear parton densities. The emphasis is on small values of x, especially 
on scale dependence of nuclear shadowing. As the main result, we show that a consistent 
picture can be obtained within the leading twist DGLAP evolution, and in particular, that 
the calculated Q 2 dependence of F^/F^ agrees very well with the recent NMC data. 
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1 Introduction 



Structure function FA has been measured in deep inelastic lepton-nucleus scatterings for 
a wide range of nuclei A [p]]. The ratios of F 2 A to the structure function of deuterium, 
F^(x,Q 2 )/F<P(x,Q 2 ), reveal clear deviations from unity. This indicates that parton dis- 
tributions of bound nucleons are different from the ones of free nucleons: xfi/A{%,Q 2 ) 7^ 
xfi/ p (x,Q 2 ). The nuclear effects in the ratio F^jF® are usually divided into the following 
regions in Bjorken x: 



• shadowing; a depletion at 2 <0.1, 

• anti-shadowing; an excess at 0.1 < x< 0.3, 

• EMC effect; a depletion at 0.3<x<0.7 

• Fermi motion; an excess towards x — > 1 and beyond. 



At the moment, there is no unique theoretical description of these effects; it is believed 
that different mechanisms are responsible for them in different kinematic regions. For an 
overview of the existing data and different models, we refer the reader to Ref. (![]. 

While the x dependence of the nuclear effects in FA was observed already in the early 
measurements Q, the Q 2 dependence is much weaker and has therefore been much more 
difficult to detect. Only recently the first observation of a Q 2 dependence of the ratio 
F2/F2 has been published by the New Muon Collaboration (NMC) ||. From the point of 
view of nuclear parton distributions, a high-precision measurement of the scale dependence 
of F2 /F® or F^JF^ is very important for pinning down the gluon distributions in nuclei 
In Ref. f| preliminary NMC data || and the small- x limit of the DGLAP equations |7j were 
used to determine the gluon ratio xg^jxgc from the measured (^-evolution of F 2 Sn /F 2 c . 

In this paper, our goal is to study whether the observed Q 2 evolution of F 2 Sn /F 2 c || is 
consistent with the leading twist DGLAP-evolution. By using the data from deeply inelastic 
lepton-nucleus scattering f8]|-||13|| and from the Drell-Yan process in pA collisions fl4 |, and 



by simultaneously requiring baryon number and momentum conservation ||, p~5|] , we will first 



determine a set of nuclear parton distributions at an initial scale Qq. When doing so, we 
want to avoid using any specific model for the nuclear effects. We then evolve the initial 
distributions up to higher scales by using lowest order DGLAP equations without parton 
fusion corrections []16 |. We will show explicitly how the nuclear effects in FA/F® and the 
similarly defined ratio for the Drell-Yan cross sections evolve in Q 2 , and make a comparison 
with the data. Our procedure leads to a consistent picture, and as the main result, we will 
show that a very good agreement is obtained between our leading twist QCD approach and 
the NMC data § for the Q 2 dependence of F 2 Sn /F 2 c . 

Scale evolution of nuclear effects has been studied already earlier |17|, [15], [§] but since 
then the parton distributions of proton have become much better known, especially at small 
values of x, where a rapid rise of the structure function F^ has been observed at HERA 



18]. Consequently, also the gluon distributions can be much better determined. So far, 



there has been no signs in the HERA data for the need of parton fusion corrections JT6[ to 
the evolution equations down to Q 2 ~ 1 GeV 2 and x ~ 10 -4 [19]. The situation has also 
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improved regarding nuclear data; more high-precision nuclear data for several mass numbers 
A are available P|-[l3| at small values of x where our main interest is focussed. The QCD 



scale evolution of nuclear effects in parton distributions, especially with specific models for 
the initial distributions at Q 2 = Ql, has been studied in O, EUl [21]]. See also Refs. |22||-|p6]. 



2 Nuclear parton distributions at Ql 

As the first task, we determine the nuclear parton distributions at an initial scale Qq as 
model-independently as possible. It turns out that some assumptions are needed about the 
nuclear effects on the initial distributions of individual parton flavours; we will try to be 
quite explicit on these assumptions. 

Fig. [I] illustrates the problem: both the deep inelastic lepton-nucleus data (DIS) and the 
proton-nucleus Drell-Yan data (DY) lie on some curves in the (x, Q 2 )-plane, determined by 
the kinematical limits and the experimental acceptances. In particular, they are not along 
a constant Q 2 line, as would be preferable for solving the DGLAP evolution. Therefore, the 
initial distributions have to be determined iteratively in such a way that the scale evolved 
distributions are consistent with the data. Like in the case of a free proton, the initial parton 
distributions at the chosen Qq serve as nonperturbative input. Our aim here is mainly to see 
whether a consistent description based on leading twist QCD scale evolution can be obtained. 
Therefore we do not attempt to implement a ^-minimization procedure in determining the 
initial distributions but plan to return to this in the future. 

The ratio of the structure function F 2 per nucleon in a nucleus A with Z protons and 
A — Z neutrons, and F® of deuterium can be written as 



F 2 A (x, Q 2 ) 
F 2 D (x,Qi 

[F 2 p/A (x, Q 2 ) + F^\x } Q 2 )] + (2Z/A - 1)[F| /A (x, Q 2 ) - F? /A (x, Q s 



F£'"(x t Q*) + F?"(x,Q 



•(2) 



In the lowest order in QCD-improved parton model F 2 (x, Q 2 ) = J2 q Cq[xq(x, Q 2 )+xq(x, Q 2 )}- 
Below the charm-mass threshold we then have 

tdA ( n 2N _ 5(ma + u A + d A + d A ) + 4sa + (jx - l)3(u A + u A - d A - d A ) , . 
K ^ x ^>- 5(u + u + d + d)+As ' (j 

where u = u(x, Q 2 ) is the known distribution of u quarks in a free proton, and u A = 
u p / A (x, Q 2 ) is the average li-quark distribution in a bound proton of a nucleus A, and similarly 
for other quark flavours. For isoscalar nuclei d n / A = u A and u n / A = d A , and in the formula 
(§) above we have assumed that this is a good approximation for the non-isoscalar nuclei as 
well. Nuclear (shadowing) effects in deuterium have been neglected in Eq. (|3]). These have 
been studied in ||28|| , where it is shown that the shadowing corrections to (F 2 + F 2 )/2 are 
of order 1 % at 0.007. Cumulative effects for x > 1 will also be neglected as well as the 
binding energy effects when using M A j A pa m p rs m n . 

We define the nuclear valence quark distributions in a usual way, q A = q A — q A , assuming 



that any differences between quarks and antiquarks in the nuclear sea 2S can be neglected. 
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Let us also define the following ratios for each sea- and valence-quark flavour: 



For later use, it is convenient to define also the corresponding ratios for the total nuclear 
valence-quark and light sea-quark distributions: 

TjAt n 2\ _ u y(x, Q 2 ) + dy{x, Q 2 ) . . 

- u v (*,Q 3 ) + dv(*,Q 2 )' w 

P A/ (o2n = "a(j, Q 2 ) + ^4^, Q 2 ) + s A (x, Q 2 ) , . 

" n(a;,Q 2 ) + d>,Q 2 )+s(x,Q 2 ) ' lJ 

and similarly for the sum w + d, and for the differences (appearing in the non-isoscalar part) 
u — d and Uy — dy 



_ u A (x,Q 2 ) 


+ d A (x,Q 2 ) 


u(x, Q 2 ) 


+ d{x,Q 2 ) 


_ u A (x,Q 2 ) 


- d A (x, Q 2 ) 


u(x, Q 2 ) 


-d(x,Q 2 ) 


u v {x,Q 2 ) 


-d v (x,Q 2 ) 


u v (x, Q 2 ) 


- d v (x, Q 2 ) 



Rt^Q 2 ) = (9) 

The relations between R v , Rt v -d v an d Rt v ->Rd v -> an d between -R^ + j, -R^_j and R?,Rf are 
obvious. Then Eq. (|3|) becomes 

>A r>2\ _ aIS(„ n2\DA/„ rfl\ 1 aISi 



R F2 (x,Q 2 ) = A^ix^Q^Ryix^Q^ + A^Q^R^^Q^+A^Q^R^Q 
,2Z 
~A 



+(— " l)[Ar(^Q 2 )Rl-dAx,Q 2 ) + AT(x,Q 2 )Rtd(^Q% (11) 



where the coefficients are known: 





%Q 2 ) 


= 5{uy(x,Q 2 ) + d v (x,Q 2 )}/N F2 (x,Q 2 ) 


(12) 


A is i 


x,Q 2 ) 


= I0[u{x, Q 2 ) + d(x, Q 2 )]/N F2 (x, Q 2 ) 


(13) 


A s [ 


x,Q 2 ) 


= As{x,Q 2 )/N F2 {x,Q 2 ) 


(14) 


A^l 


>,Q 2 ) 


= 3[u v {x,Q 2 )-d v {x,Q 2 )]/N F2 {x,Q 2 ) 


(15) 


aNIS, 
n ud 1 


>,Q 2 ) 


= 6[u(x, Q 2 ) - d(x, Q 2 )}/N F2 (x, Q 2 ) 


(16) 


N F2 { 


x,Q 2 ) 


= 5[u v (x, Q 2 ) + dvix, Q 2 )] + 10[u(x, Q 2 ) + d(x, Q 2 )] + As{x, Q 2 ). 


(17) 



In order to have further constraints for the individual nuclear ratios above, we will con- 
sider differential Drell-Yan cross sections in pA collisions. For a nucleus with Z protons and 
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A — Z neutrons the ratio of cross section to that for deuterium can be written in the lowest 
order asf] 



, | , ,-,2 ji^C^Y I dx 2 dQ 2 



dy 2,Q ^da^y I dx 2 dQ 2 

= {A[ui(u A + d A ) + u x {u A + d A )} + [dtidi + u A ) + d x {d A + u A )} + 4 Sl s A + ...}/N DY 
17 

+(— - l){A[ Ul {u A - d A ) + u x {u A - d A )] + [di(d} - u A ) + d x {d A - vim/Nor (18) 
where 

N DY = 4[mi(m 2 + d 2 ) + ui{u 2 + d 2 )\ + [di{d 2 + u 2 ) + di(d 2 + u 2 )} + 4sis 2 + ... (19) 

and we have used the notation q\ A ^ = Q 2 ) for i = 1,2 and q = u,d, s, .... The variable 

Q 2 is the invariant mass of the lepton pair. The target (projectile) momentum fraction is 
x 2 (xi), and X\ = Q 2 /(sx 2 ). With the definitions in Eqs. (|4T)-(|i~0|), we obtain 

Ri Y (x,Q 2 ) = Bi s d (x 1 ,x 2 ,Q 2 )R A +d (x 2 ,Q 2 ) + B v s (x 1 ,x 2 ,Q 2 )R A (x 2 ,Q 2 ) + 



?/ 

-B s ( Xl ,x 2 ,Q 2 )R A (x 2 ,Q 2 ) + (— - l)[B^ d IS (x 1 ,x 2 ,Q 2 )R A _ (1 (x 2 ,Q 2 ) + 



,2Z 
' A 

+B$ IS (x 1 ,x 2 ,Q 2 )R A v _ dv (x 2 ,Q 2 )} (20) 
where the coefficients are: 



B'Jix^x^Q 2 ) = [A(u 1 +u 1 ) + d 1 + d 1 }(u 2 + d 2 )/N DY (21) 

Bl?(x 1} x 2 ,Q 2 ) = (4u 1 + d 1 )(u V2 + dv 2 )/N DY (22) 

B s (x l} x 2 ,Q 2 ) = A Sl s 2 /N DY (23) 

B% d IS ( Xl ,x 2 ,Q 2 ) = [A(u 1 +u 1 )-d 1 -d 1 ](u 2 -d 2 )/N DY (24) 

B$ IS ( Xl ,x 2 ,Q 2 ) = (Au 1 -d 1 )(u V2 -d V2 )/N DY (25) 

N DY = [A(ui + ui) + di + d~i]{u 2 + d 2 ) + Asis 2 + (Aui + di)(u V2 + d V2 \2&) 

with q V2 = qv{x 2l Q 2 ) and q x = qy(xi, Q 2 ) + q(xi, Q 2 ), q = u,d. 



The DIS data ||, |TU], [XT], [3] which we use in determining the nuclear ratios, are ap- 
proximately corrected for non-isoscalar effects, so we may consider isoscalar nuclei first. In 
Eqs. ( |Tl"D and ( |20"D the terms proportional to (2Z/A — 1) can be dropped, and we are left 
with Ry, R A +d and R A only. Without additional information we cannot fix three ratios from 
two equations. Further constraints could in principle be obtained from the measurements 
of nuclear structure functions F^/F^ and F^ A /F^ D with neutrino beams [30], but more 



statistics and mass number systematics would be needed. Instead, as a first approximation, 
we take for the sea quarks R A (x, Ql) = R A +d (x, Ql). For evolving each flavour separately, we 
will also assume that initially R A (x,Q$) = R A (x,Ql), and similarly for the valence quarks 
R a v (x,Qq) = R^^x^Qq) = R A (x,Ql). It should be emphasized that this approximation is 



4 In Rcf. [Q, the corresponding formula is incorrect, but it has the correct large-x_F limit, which was used 
only. 
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needed only at the initial scale Ql, in determining the initial distributions for the DGLAP 



evolution. In this paper we use, for simplicity, the parton distribution set GRV-LO [[31 
where the massive quarks evolve as massless ones above each mass-threshold. This set also 
assumes u = d, so R£ = Rf at all scales for this particular set. 

In this approximation, we obtain two equations for the ratios at an initial scale Ql, 

Ri(x,Q 2 ) = A^{x,Ql)R^{x,Ql) + [A^(x,Ql)+A s {x,Q 2 )]R^(x,Ql), (27) 
Rdy&QI) = B^ix^x^^R^Ql) + [Bil( Xl ,x,Ql) + B s ( Xl ,x,Ql)]Rj(x,Ql) (28) 

where Ql is chosen below the charm threshold. Eqs. (|27j ) and (|28| ) would fix R v and R^ if 
the DIS and DY data would lie at the same values of x and Ql- Since this is not the case 
(see Fig. [I]), the initial profiles R$(x,Ql) and R v (x, Ql) are determined iteratively after 
comparing the evolved distributions with the data. As is clear from Eqs. (^) and (28), the 



initial nuclear ratios acquire some dependence on the particular parton distribution set used 
for the proton. With the modern sets, however, this dependence can be expected to be quite 
weak. 

As a further constraint for nuclear valence quarks ffl, baryon number conservation 



f 1 f 1 

/ dx[u v (x, Ql) + d v (x, Ql)]Ry(x, Ql) = / dx[u v (x, Ql) + d v (x, Ql)] = 3. (29) 

JQ J 

will also be required |1| §]. 

In practice, our iteration procedure for determining the ratios R$(x, Ql) and R v (x,Ql) 
proceeds with the following steps: 

• make an ansatz for Rp 2 (x, Qq) based on the DIS data; 

• decompose Rp 2 (x, Ql) into Ry and R$, and constrain R v with baryon number conser- 
vation; 

• estimate Rq(x,QI) (see the discussion below); 

• perform DGLAP evolution with the obtained initial nuclear parton distributions; 

• constrain R$(x, Ql) and R v (x,Ql) with the DY data. 

First, we introduce an initial parametrization for Rp 2 (x, Ql). We choose Ql = 2.25 GeV 2 , 
which is the charm-mass threshold for the GRV-LO set. Rather than trying to make a 
separate analysis for each nucleus, we parameterize also the A dependence of Rp 2 . The 
functional form is given in the Appendix. Note also, that as the existing parametrizations 
p2| , [33j are fits to the data along the kinematical curves in Fig. [I], we cannot use them 
directly for obtaining the distributions at fixed scale Qq. Furthermore, we have to make 
an assumption on the behaviour of Rp 2 (x, Ql) at x< 10~ 3 . In this region a saturation of 
shadowing has been observed [12|, || at nonperturbative scales, (Q 2 ) Ql (see Figs. [I] and 



0). Motivated by this, we will also assume saturation of shadowing in Rp 2 (x, Ql), keeping 
in mind, however, that consistency of this assumption must be verified after evolving the 
distributions (see Figs. £|). 
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We then decompose Rp 2 into Ry and R$ according to Eq. fl2"7|). First observation is 
that we cannot have R v = Rp 2 at all values of x, because up to 9 % of baryon number 
would be missing for the heaviest nuclei. Second observation is that Rp ~ R v (R$) at large 
(small) values of x. In practice, in the region of large (small) x, it becomes impossible to 
determine R$ (Ry) directly from the data for Rp 2 , because of the negligible sea (valence)- 
quark contribution. Therefore, we apply a piecewize construction: at x < x p ~ 0.1 we fix 
Ry with the same functional form as for Rp 2 but with different parameters. Then, through 
Eq. ( p?D Rg becomes fixed. At x p < x < x eq ~ 0.4 we fix R$ in turn, with a simple form of 
a plateau Rs(x p < x < x eq , Ql) = Rg(x p , Ql). Now the approximate plateau in R$ controls 
the height of the anti-shadowing peak in R v . At x = x eq the sea-quark ratio Rg(x p ,Ql) 
becomes equal to Rp 2 (x p , Qq), and at x > x eq , in lack of further information on Rg, we again 
use the simplest approximation Ry = Rg = Rp 2 . The precise values of x p and x eq , (i.e. 
the location and height of the plateau in R$), together with the parameters for Ry are first 
constrained by the baryon number conservation at Q 2 = Ql and after the DGLAP evolution 
by the DY data 0. 

Finally, we define the nuclear gluon ratio by 

R^x,Q 2 ) = g A (x,Q 2 )/g(x,Q 2 ), (30) 

and specify Rq(x, Ql) to have the full input for the DGLAP evolution. We constrain 
Rq(x, Ql) with momentum conservation [jT5|, p|: 

1 = J dxx{g{x, Ql)R G (x, Ql) + [u v {x, Ql) + d v {x, Q 2 )]R${x, Ql) + 

2[u(x, Ql) + d(x, Ql) + s(x, Q 2 )]Rj(x, Q 2 )|, (31) 

where Ry and R§ are determined as described above. Compared to free nucleons, we find 
that some momentum is transferred from quarks to gluons in nuclei. This effect is not very 
large: in A = 208 the glue carries about 4 % more momentum than in a free nucleon. This 
is in agreement with the earlier studies [13, ffl. Consequently, without any x-dependence in 
Rq(x,QI) we would have Rq(x,QI) = 1.04. 

Since the sea quarks are shadowed at small x, we expect shadowing of the nuclear gluons 
as well. A requirement of stable scale evolution can be used together with the recent NMC 
data to further constrain nuclear gluon shadowing. At the small-x limit of the DGLAP 
equations one obtains 

dR^Q 2 ) dF 2 D (x,Q 2 )/d\ogQ 2 f dF 2 A (x,Q 2 )/dlogQ 2 A 2 1 
dlogQ 2 ~ F 2 D (x,Q 2 ) \dF 2 D (x,Q 2 )/d\ogQ 2 WVJ J 1 ' 



5a s xg(2x, Q 2 ) 
~9^F 2 D (x,Q 2 ) 



R£(2x,Q 2 )-Rp 2 (x,Q 2 )y (33) 



where we have used the result dF 2 (x, Q 2 )/d\ogQ 2 ~ 5a s xg(2x, Q 2 )/9ti, derived in P^j . As 
for the initial Rp 2 (x,Ql), we will also assume a saturation in shadowing of the glue, so 
that R£(2x,Ql) w R£{x,Ql) at the limit x -»• 0. For the GRV-LO set we are using, the 
ratio xg(2x, Ql)/F 2 (x, Ql) is non-zero and grows slowly at small values of x, so Rp 2 evolves 
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towards Rq due to the factor Rq — Rp 2 in Eq. (|33|). An initial condition stable in the 
evolution is obtained by requiring Rq(x, Qq) ~ Rp 2 {. x i Ql) a ^ ver y sma U values of x. 

Due to momentum conservation, the loss of gluons in the shadowing region must lead 
to an increase of gluon distribution somewhere at larger x. In Ref. [|J, it was assumed that 
Rq(x,Q 2 ) is constant at large x. This, however, resulted in somewhat unstable profiles for 
Rq and R$ in the evolution. Therefore, we expect that there is an EMC-effect also for the 
glue, and that Rq(x, Ql) is peaked around the same values as Ry(x, Ql) and R$(x, Ql). 

In practice, we determine the gluon ratio as follows. We start with Rq(x, Ql) ~ Rp 2 (x, Ql) 
at small values of x. Should we use this equality for all values of x, some momentum would 
be missing for all nuclei (11 % for A = 208). This indicates that conservation of momentum 
requires quite strong ant i- shadowing for the gluons. For Rq(x, Ql) we use the functional form 
given in the Appendix. In determination of the position and width of the anti-shadowing 
peak (which we assume to be independent of A) we use Ref. to constrain the value of x 
where xgsn/xgc ~ 1- The amount of anti-shadowing then follows from momentum conser- 
vation, leading to a good agreement with Ref. || once the position and width parameters of 
the peak in R G are fixed. 

The initial nuclear ratios R G (x,Ql), R v (x,Ql) and R$(x, Ql) at Ql = 2.25 GeV 2 are 
shown for isoscalar nuclei in Fig. ||] together with the ratios Rp 2 (x, Ql). Correlations in these 
ratios are easily understood from the figures: the more R v is shadowed, the more it has 
anti-shadowing, and the more R$ is suppressed in the anti-shadowing region. Similarly, due 
to momentum conservation, the more shadowing Rq has, the more anti-shadowing is needed. 
As seen in the figures, there is more anti-shadowing and less shadowing for Ry than for Rp . 
Correspondingly, R$ is more shadowed and not anti-shadowed at all at Qq = 2.25 GeV 2 . The 
gluon ratios R G turn out to be clearly more anti-shadowed than Rp 2 and Ry. In comparison 
with the earlier studies, the lack of sea-quark anti-shadowing agrees with JL5], £§]. We get, 



however, somewhat more anti-shadowing for gluons than in []T3, which is a consequence 
of the small- X enhancement of gluon densities in a proton. 



3 Scale evolution and results 
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arc 



Scale evolution is straightforward to carry out once the parton fusion corrections 
neglected. In the HERA data for F% (x, Q 2 ) [|18|, [19] there is no evidence of the fusion cor- 
rections at Q 2 <; 1 GeV 2 and x ^ 10~ 4 . For nuclei, parton fusion should be stronger due to its 



expected A 1 / 3 scaling |IBL E5[, but since we start the evolution at Ql = 2.25 GeV 2 > 1 GeV 2 



the effects of these corrections should be small 0], at least in the x range of the NMC data 
0j. The role of the fusion corrections is, however, a very interesting question [ff7|, f|, ^0, ^3 



which deserves further analyses with constraints for the proton from the new HERA data 
at very small Q 2 and very small x. 



We use the parametrization of GRV-LO [31 1 for the initial parton distributions of free 



protons. We evolve these and the initial nuclear distributions from Qq = 2.25 GeV 2 to 
Q 2 ~ 10000 GeV 2 by using the DGLAP equations for gluons and valence and sea quarks of 
each flavor. Massive quarks are generated through the evolution above each mass threshold 
as in |^TJ. We solve the evolution equations directly in (x, Q 2 )-space. The resulting scale 
evolution of the ratios Rq, R$ and Ry is shown in Figs. |3| for A = 208. Compared with an 
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earlier study |4[] , the difference in the evolution of Rq and Rp 2 is not as dramatic, due to the 
enhanced, more stable gluon and sea-quark distributions at small vales of x. 

In Figs. f| we plot the scale evolution of the ratio Rp 2 (x, Q 2 ) for isoscalars ^He, 1 gC, ^Ca. 
and for an artificial isoscalar nucleus A = 208 at different fixed Q 2 . The data shown and 
used in the analysis, are the re-analyzed data from NMC g and SLAC 0. For a more 
transparent comparison, we show our calculation with filled circles at Q 2 equal to the (Q 2 ) 
of the NMC data at different x. 

In Fig. |5|, we show the behaviour of Rp 2 (x, Q 2 ) at very small x, together with the NMC 
data || and || for 12 C. This figure confirms the consistency of our initial assumption of 
the saturation of shadowing in Rp 2 (x, Qq): at fixed x, shadowing decreases with increasing 
Q 2 and therefore our initial Rp 2 (x, Qq) should not lie below the data measured at Q 2 < 
Qq. This of course implies an implicit assumption that the scale dependence of Rp in 
the nonperturbative region is towards the same direction as in the perturbative region. It 
should be kept in mind, however, that at Q 2 > Qq the sign of dRp 2 (x, Q 2 )/dQ 2 , as indicated 
by Eq. (|33|) , actually depends on the gluon shadowing at small x, i.e. if the gluons were 
more shadowed, the decrease of shadowing in Rp would be slower, or shadowing could even 
increase. In this case Rp 2 (x, Ql) should be re-determined accordingly. In this way, as pointed 
out in [||, at small values of x the scale evolution of Rp 2 reflects the amount of nuclear gluon 
shadowing. 

Small- x data for Rp exist also from the E665 collaboration [13], but they lie above the 



NMC data (see e.g. |10[| ). Because of the smaller error bars in the NMC data, we decided not 
to use the absolute E665 data for determining Rp 2 (x,Ql) here. However, as pointed out in 



TP ], if one considers the ratios of different nuclei, the two data sets are consistent. Therefore, 
we have also used the ratio F^ h / -F 2 as obtained from the E665 data [13) . In Fig. |6] we show 
the results at small x as in Fig. |5| but for the ratio of ratios, Rp 2 / Rp 2 with the data from |10 



and [|HJ. Note that the inner (outer) error bars under(over)estimate the real errors of the 
E665 data (see the figure caption). Even without a more precise x 2 minimization, this figure 
partly helps us in constraining the A dependence at small values of x in the parametrization 
of Rp 2 given in the Appendix. 

The systematics in nuclear mass number A, especially at x<0.l, is presented by the 
NMC collaboration in and we have also made use of these data. In Figs. [7] we show the 
scale evolution of F^/F^, i.e. of Rp 2 /Rp 2 , and the comparison with the NMC data. Since 
the data are (approximatively) corrected for non-isoscalar effects, we have always A = 2Z. In 
our analysis, the data for the (A = 117)/C ratio gives the most stringent constraints. Again, 
for an easier comparison, the filled symbols show our calculation at the (Q 2 ) of the data at 
different x. 

How the Drell-Yan data from the E772-collaboration [H enters our analysis, can be seen 
from Figs. |8| where we have plotted the ratio of Drell-Yan cross sections as defined in Eq. 
( POD at different invariant masses Q 2 for ^C, 2oCa, ggFe an d ^IW. The non-isoscalar terms 
for Fe and W have been taken into account according to Eq. (|20"D, by using the corrections 
obtained with isoscalar nuclei. As before, the filled symbols show our calculation at the mass 
values (Q 2 ) of the data J27J. The constraint for determining the initial ratios Rg(x, Qq) and 
Ry(x, Ql) enters as follows: let us suppose that valence quarks are shadowed less than in 
Figs. This leads to a smaller anti-shadowing for Ry(x, Ql) because of baryon number 
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conservation, Eq. fl2"9"|). At the same time, Eq. (|3]) brings -Rg(x, Ql) closer to Rp at x ~ 0.1, 
i.e. Rg becomes less suppressed. In the scale evolution, there is a cross-over region in Rp, Y 
where evolution is very weak (see Figs. and the location of which depends on the plateau 
in Rg(x ~ 0.1, Qq). With less suppressed Rg at x ~ 0.1, the location of the cross-over region 
shifts to larger x. Consequently, the evolution of Rp, Y becomes too fast at small values of 
x, and the calculated values overshoot the data. On the other hand, if R v is much more 
shadowed than in Figs. [2], the cross-over region in R^ Y shifts to the left, scale evolution 
becomes too fast at x £ 0.1, and the calculation falls below the data. Iterating between these 
possibilities we have obtained the initial R$(x, Qq) shown in the figures. 

In Fig. |9|, comparison with the results of Ref. || for the gluon ratio gsn/ 9c is shown. 
Again, the curves are plotted at fixed values of Q 2 and the filled symbols show our calculation 
at the (Q 2 ) of the preliminary NMC data || as quoted in |5[]. It should be noted that even 
though only the point where gsn/ gc ~ 1 was constrained to agree with ||, the resulting 
amount of shadowing and anti-shadowing comes out in a very good agreement with the 
results of Gousset and Pirner. 

Finally, in Figs. 0, we present the main result of our study, comparison of the calcu- 
lated scale evolution of F 2 /F 2 with the recent NMC data ||. The data is plotted with 
statistical errors only because we are more interested in the slopes of Q 2 evolution than in 
the absolute normalization to the data. Notice that normalization of our curves at each 
fixed x in Figs. [10] depends on the initial ratio shown for Sn/C in Fig. [7[ In fact we obtain 
a fairly good agreement also for the overall normalization of the data. The changes in the 
sign of <9(_F 2 Sn / 'F^) I dQ 2 do not seem to be in contradiction with the data, either. It is also 
interesting to notice that at small values of x the slope is not linear in logQ 2 . 

After performing the scale evolution, let us look back into the assumptions made in 
constraining the nuclear ratios at the initial scale Q^. As expected on the basis of Eq. fl33|) , 
the approximation Rq = Rp 2 for the initial distributions at extremely small x turns out 
to be quite stable in the evolution; at the smallest x we consider, x = 10~ 6 , deviations 
from R£ = R F2 are only about 5 % up to Q 2 = 10000 GeV 2 for A = 208. For the valence 
quarks, the initial approximations R^ v = R v (x,Ql) and R^ y = Ry(x,Ql) are very stable: 
deviations are within 1 % at x < 0.7 at all scales considered. For the GRV-LO set IBII, there 



is no difference between u and d, so once the initial approximation R^ = R^ is made, it 
holds exactly at all scales. The initial approximation R^ + ^ = R$ is also a stable one, the 
deviations are within 2 % during the evolution. Initial approximation Rf = R$ is slightly 
less stable, but deviations do stay within 7 % at all scales considered. The deviations are 
concentrated in the region x = 0.01. ..0.1, where the nuclear effects for sea quarks are small 
in any case. During the evolution, the strange quarks tend to develop more anti-shadowing 
than Rg, while the u and d have an opposite tendency relative to R$. More work is needed 
to relax the assumption R^(x,Ql) = R^ + j(x,Qq) but this additional correction is beyond 
the scope of this study. 

4 Discussion and Conclusions 

The main emphasis of this study is in the (Revolution of nuclear effects in parton distri- 
butions in the region of small x. We have used the deep inelastic scattering data from IA 
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collisions [p||-||13|| and the Drell-Yan dilepton data from pA collisions |T3| together with con- 
servation of baryon number and momentum to constrain the initial nuclear distributions as 
model-independently as possible. As the main result, we have shown that a consistent pic- 
ture arises, and a very good — almost surprisingly good — agreement with the measured Q 2 
evolution of the structure function ratio _F 2 Sn /.F 2 c |§ can be obtained already with the lowest 
order leading twist DGLAP evolution. We can therefore conclude that the effect of parton 
fusion corrections in the evolution |1| is neg ligible, at least at x > 0.01 and Q 2 > 2.25 GeV 2 . 
In this region, the nuclear modifications are effectively built in the nonperturbative initial 
conditions for the DGLAP evolution. This result also agrees with Ref. |p0[ . 

We point out, however, that even though we obtain a very good agreement with the 
NMC data || and with the analysis of Ref. ||, we can confirm our initial assumption of 
gluon shadowing at small values of x only on fairly qualitative grounds (stability of the 
evolution), rather than through a direct comparison with the data. The reason for this 
is seen qualitatively from Eq. (|33|) : the data on F 2 Sn /F 2 c at x = 0.0125 constrains R G at 
x = 0.025, which is only the beginning of the gluon shadowing region and where nuclear 
effects in Rq are not very strong (see Figs. |2|). To constrain the gluon shadowing further, it 
would be very important to have data for the Q 2 dependence of F 2 Sn / F 2 available at smaller 
values of x. 

The NMC data on deep inelastic J/\l/-production in Sn/C |36[ is not included in our 
analysis but it is interesting to notice, as pointed out in ||, that the amount of gluon 
ant i- shadowing in the ratio Sn/C agrees with the excess reported in |36|j . The gluon anti- 
shadowing we obtain is also consistent with the E789 data on D meson production in pA 
collisions [37], although the error bar on the data point is quite large. With the strongly 
interacting final states, higher twist production mechanisms |38[ may well make the picture 
more complicated regarding factorization. However, if the initial state effects can be factor- 
ized into the nuclear parton densities [39|, our analysis on R G should be of direct use also 
for computing strongly interacting final states in nuclear collisions. Our results should also 
provide more insight in understanding the recently measured increase in J/\l/ suppression 
relative to the Drell-Yan background in central Pb-Pb collisions at the CERN SPS f40 |. 

Nuclear structure functions F% Fe ^ F^ Fe Ell] are used in the global analyses of parton 



distributions for a free proton [ |2| . Our results should offer a more consistent way of unfolding 
the nuclear effects there.f] Also the coordinate space description of nuclear parton densities 



43] could be studied in more detail by using the scale evolved nuclear ratios we have presented 



here. Related to the coordinate space description and scale evolution of nuclear parton 
distributions, the connection of our results with those in the region of very small x and very 



large A where the higher twist terms can be expected to be more important |[26| , should also 
be studied in more detail. 

To conclude, the results of our study are encouraging. They show that further QCD- 
analysis on nuclear parton distributions is worth performing. For a more detailed analysis, 
this study can be improved in obvious ways: A more quantitative error analysis should be 
done by performing a minimization of x 2 i n fits to the data. In determining the nuclear 
gluon distributions by using the NMC data || in a more consistent manner such analysis 
would be required. Also more work remains to be done in determining the nuclear ratios 



5 We thank C.A. Salgado for discussions on this point. 
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for individual quark flavors. This should be done, however, model-independently by using 
the measured data wherever possible. More systematics on mass number dependence of 
the structure functions F£ A and F£ A and increased statistics against deuterium would be 



helpful in extracting the nuclear ratios of individual parton distributions (see also [fH]]). Also 
data on F A /F® (or e.g. F 2 A /F^), not corrected for non-isoscalar effects, might be useful 
for constraining the difference between R A and Rf, provided that precision of the data is 
sufficient. Naturally, further measurements of the Drell-Yan cross sections in pA collisions 
with more statistics would be very useful. 

We do not expect the nuclear ratios to depend strongly on the choice for the (modern) 
parton distributions of the free proton, but an explicit study of this is in progress. For 
numerical applications we are also preparing a package which produces scale dependent 
nuclear ratios Rf(x, Q 2 ) for any parton flavour / in an arbitrary nucleus A |[45|| . 

Eventually, our analysis should be extended to next-to-leading order in the cross sections 
and in the scale evolution^]. It will also be interesting to study the role of the parton fusion 
corrections in more detail in the light of the HERA results. 
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Smirnov, M. Strikman, R. Venugopalan and M. Vanttinen for useful discussions, and P. 
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during the Hard Probe meeting in September 1997, where part of this work was discussed 
and developed. We also thank the members of the Hard Probe Collaboration for helpful 
discussions, comments and encouragement. This work was supported by the Academy of 
Finland, grant no. 27574. 

Appendix 

We use the following piecewize parametrization for the ratio Rp 2 (x, Qq), motivated by 



R F2 (x,A) = 1 + s(x,A) +emc(x,A) + f(x,A), (34) 
with the shadowing part given by 

s(x, A) = [s (x, A) - l]e- x2 / x °Q{x m - x), (35) 

where 6 is a step function, and 

/ a \ _ 1 +a s k 2 (~L/x- l/x se )p(x) 
* J "l + «^(l/x-l/x se M,)' [6b) 

and p(x) = max[l, (x/x se ) m } and p$ > 0. 
The emc-part is controlled by a function 

emc(x, A) = (a + a x x + a 2 x 2 + a 3 x 3 - 1)(1 - e~ {x/x ^ )2 )Q(x m - x), (37) 



3 NL0 evolution of nuclear parton distributions is studied in [[20) 
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Table 1: The parameters used for Rp (x,Qq) and R v (x < x p ,Qq) with x p = 0.09 (see Sec. 2). 



where x$ e = XQ e (A/12) Pe . We fix the location of the minimum in Rp 2 (x, A) at x m = 0.74 for 
all nuclei. With the parameter x for s(x, A) and with the parameter Xq £ for emc(x, A), we 
obtain a smoothly behaving anti-shadowing region in Rp 2 (x, A). 
The Fermi-motion region is parametrized with 



f(x,A) 



'1 _ jA\ \ X X m) j. 

k J ml i \2 'J i 

\Xf X m J 



G(x 



(3* 



with = /^ a (A/40) a , where a is given by Eq. (9) of Ref. [|TTJ with x = x m . For simplicity, 
we have fixed Rp 2 (xf, A) = 1 for all nuclei. 

The parameters ao, ax, 02, 03 above are determined from the following conditions: 



e(x e ) 



e(x m ) = /, 



,4 

m ' 



<9e(x) 



<9x 



0. 



<9e(x) 



x=x 0e 



dx 



0. 



(39) 



where e(x) = + aix + CI22 2 + c^x 3 . In the Table [I] below, we give the parameters for the 
parametrizations of Rp 2 (x, Ql) and R v (x < x p ,Ql). 

The form of the initial nuclear gluon ratio Rq(x, Qq) is the following: 



R£(x,Q 2 ) = R F2 (x,A)\l + N 



exp 



(log a; - (\ogx g - Ap/1.5)) : 



+ 



exp 



(logs - (logXg + Ag/1.5)) 

^9 



2\ n 



(40) 



where x 9 gives the position for the additional anti-shadowing bump in Rq, and A g determines 
its width in log x. The amplitude N for the enhancement can be directly determined by using 
the momentum sum rule (|3T|). In Figs. |2|, we have x g = 0.09 and A g = 0.9. 
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Figure 1: Typical correlation of the scale (Q 2 ) and x in measurements of F^(x,Q 2 ) in deeply 
inelastic I A scatterings and correlation of the invariant mass (Q 2 ) and x = X2 of Drell-Yan cross 
sections measured in pA collisions. The correlations in some of the NMC data |8| (solid lines), [fLOf 
(dotted-dashed), Q (dotted) and in some of the E665 data [12] (dashed), and in the E772 data 
|l4], [27j (long dashed) are shown. The horizontal dotted line illustrates the initial scale Qq we have 
chosen and above which we perform the DGLAP evolution of nuclear parton densities. 
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Figure 2: The initial nuclear ratios Rq(x,Qq) (solid line), Ry(x, Qq) (dotted) and R^(x,Q%) 
(dashed) for isoscalar nuclei at Qq = 2.25 GeV 2 . The ratio Rp 2 (x, Qq) (dotted-dashed) is also 
shown. 
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Figure 3: Scale evolution of the ratios R^(x,Q 2 ), R$(x,Q 2 ), R v (x,Q 2 ) and Rp 2 (x,Q 2 ) for an 

isoscalar nucleus ^4=208. The ratios are shown as functions of x at fixed values of Q 2 = 2.25 GeV 2 

(solid lines), 5.39 GeV 2 (dotted), 14.7 GeV 2 (dashed), 39.9 GeV 2 (dotted-dashed), 108 GeV 2 

(double-dashed), equidistant in logQ 2 , and 10000 GeV 2 (dashed). For R 
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Figure 4: Scale evolution of the ratio Rp 2 (x,Q 2 ) for isoscalar nuclei A=4, 12 and 40. As in 
Fig. g, the ratios are plotted as functions of x but with Q 2 fixed to 2.25, 3.70, 6.93, 12.9, 24.2 GeV 2 , 
equidistant in logQ 2 , and 10000 GeV 2 . The reanalyzed NMC data || is shown by the boxes, the 
reanalyzed SL AC data by the triangles [|lT] . The statistical and systematic errors have been added 
in quadrature. The filled circles show our calculation at the (Q 2 ) values of the NMC data. Notice 
that the vertical scale of each panel is different. 

19 




Figure 5: The ratio Rp 2 (x,Q 2 ) shown together with the renanalyzed NMC data || (boxes) and 
the combined NMC data |J (triangles). The calculated Rp 2 are shown at the same fixed values of 
Q 2 as in Fig. Notice that at x < 0.01 the scales (Q 2 ) of the data are less than our Qq = 2.25 GeV 2 . 
The statistical and systematic errors of the data are added in quadrature. 
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Figure 6: The ratio R^(x, Q 2 )/Rp 2 (x, Q 2 ) as a function of x at the same fixed values of Q 2 as in 
Fig. The NMC data for Pb/C (T^j is shown by the squares and with statistical and systematic 
errors added in quadrature. The ratio of the E665 data for Rp 2 and Rp 2 jnj is shown with the 
triangles. The inner error bars are obtained by including the independent statistical errors only, 
and the outer ones by adding first the statistical and systematic errors separately for Pb/D and 
C/D in quadrature, and then taking these errors to be independent. This figure shows how our 
calculation relates to the small- x region where the measurements are at nonperturbative scales. 
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Figure 7: The ratios Rp 2 (x,Q 2 )/Rp 2 (x,Q 2 ) as functions of x for isoscalar nuclei. The values of 
Q 2 are the same as in Fig. ||. For comparison with the NMC data |l(| (boxes), the filled circles 
show our calculation at the {Q 2 ) of the data. The inner error bars stand for the statistical errors 
only, the outer ones for statistical and systematic errors added in quadrature. 
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Figure 8: The ratios of differential Drell-Yan cross sections in pA and pD as functions of x = X2 
for l $C/D, fjCa/A %Fe/D and ^fW/D. Our calculation for R^ Y (x,Q 2 ) of Eq. @ is shown at 
fixed values of the invariant mass Q 2 = 2.25 GeV 2 (solid line), 24.2 GeV 2 (dashed), and 139 GeV 2 
(dotted-dashed). The data shown by the boxes is from E772 [14]. In the graph the statistical 
errors and the quoted 2 % systematic errors are added in quadrature. The filled circles show our 
calculation at the (Q 2 ) of the data 
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Figure 9: The gluon ratios Rq l (x,Q 2 )/Rq(x,Q 2 ) as functions of x at fixed values of Q 2 = 2.25, 
3.27, 5.39, 8.89, 14.7 GeV 2 , equidistant in logQ 2 , and 10000 GeV 2 . The filled circles show the 
comparison with the results of Ref. ]5J , presented by the boxes. 
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Figure 10: The calculated scale evolution of F^ix, Q^/F^fx, Q 2 ) compared with the NMC data 
|P] at different fixed values of x. The data are plotted with statistical errors only. 
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